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We introduce a basic model for human mobility that accounts for the different dynamics arising 
from individuals embarking on short trips (and returning to their home locations) and individuals 
relocating to a new home. The two modes of motion occur on widely separated time scales and affect 
different phenomena: for example, short trips, constituting the bulk of human motion, dominate the 
spread of diseases, while the spread of genes or family names is more strongly affected by relocations. 
In many everyday life situations the two modes of motion need to be considered simultaneously. 
This is illustrated by a simple SIR model for epidemic spreading where trips and relocations lead 
to dramatic differences even under equal volumes of travel. 
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Wc introduce a basic model for human mobility that 
accounts for the different dynamics arising from individu- 
als embarking on short trips (and returning to their home 
locations) and individuals relocating to a new home. The 
differences between the two modes of motion comes to 
light on contrasting two recent studies, one tracking the 
geographical location of dollar bills [1], the other that of 
mobile cell phones [2] . Trips introduce two characteristic 
time scales; the time between trips, 9, and the duration 
of each trip, r, and relocations introduces a third time 
scale, T, for the time between relocations. In practice, 
T ^ years, 9 ~ months, and r days, so the three time 
scales are widely separated. Traditionally, studies incor- 
porating human motion assume only a single mode, using 
a generic rate to account for all types of motion. 

In what follows, we first introduce a social mobility 
model, capable of incorporating both modes of motion 
— trips and relocations — and show how it can be cast 
as a Markovian process, thus simplifying its analysis. We 
then argue that the distinction between the two modes is 
more that simple semantics or a minor increase in resolu- 
tion, but may have significant consequences for dynamics 
supported by the motion, such as the spreading of an epi- 
demic. We illustrate though a specific example where an 
infectious disease, evolving according to the susceptible- 
infected-recovered (SIR) model, and spreading between 
two populations through trips and relocations, exhibits 
dramatic differences in the endemic levels of disease, de- 
pending on the mode of travel. 

We now describe our social mobility model (SMM) , for- 
mulated as a Markov process. Assume that the human 
population may be viewed as a number of homogeneous 
agents, M , moving independently among N discrete ge- 
ographic locations. Because agents move independently, 
we need only describe the process for a single individual. 
To support preferential return to home, we define states 
of the system as ordered pairs 



Denote the fraction of individuals at state («, j) by T^ijit) 
and require, accordingly. 



(2) 



Under this conceptual description. 



j i 

arc the fraction of people actually at location i, and the 
fraction of people whose home is j, respectively. Note 
that TT^j would be expected to match census data, aj, 
taken at some time close to the start of the process (at 
time t = 0). We may in fact assume that the census is a 
fair representation of the equilibrium distribution, and 

aj « 7r*j(0) « n^jit), 

where the approximations are based on the idea that for 
a large population, the stochastic deviation from equi- 
librium should be small. Implicit is also the assumption 
that TTij <C (Tj for all i ^ j. 

The SMM allows for two types of motion; relocation to 
a new home and short trips. Relocations are relatively 
rare, and we denote the typical time between any two 
relocations (for a given agent) by T. Trips are more fre- 
quent and are characterized by the time between trips, 9, 
and the typical time t spent in one location while travel- 
ing. In practice, these three time scale are well separated, 
T ^ 6* ^ r, as T years, 9 ^ months, and r ^ days, 
though this is not required by our model. 

An individual at his home location, i, [state (i, i),] may 
decide to relocate to j, [state («,j),] at rate LOij/T. The 
factor LOij expresses the likelihood to relocate from i to 
j (as opposed to some other location) and is normalized 



1, such that l/T is approximately the 



= (current location, home location) 
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total rate for relocations. Rather than relocating, he may 
decide to take a trip to location j, — )■ at rate 

h'ij/9, where the i^ij indicate the preference to travel from 
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i to j and are normalized as J^i J2j '^i'^ij = li such that 
the total rate for trips is approximately 1/6. 



A person away from home, j, at state («, j), would stay 
at i an average length of time t, whereupon he might re- 
turn home, to with probability r (and rate r/r), or 
continue on his trip, to {k,j), with probability 1 — r (and 
rate (1 — r)vi)^lT). Note that in this fashion r controls 
the length of a trip: if r = 1 all trips are merely to a 
destination and back (no layovers), but more generally 
the average number of legs in a trip is (1 + r)/r. 



Finally, a person might decide to relocate while away 
from home, (z, j) — ?■ (i, fc), at rate ujj]^/T, but since tt^ j- ^ 
CTj this is an uncommon event and not much is lost by 
neglecting this possibility. We list it as an optional move. 
The various transitions allowed by our model, and their 
rates, are summarized in Table I. 



The Wy's and f^-'s are assumed as input, though we 
acknowledge that obtaining specific values is an art in its 
own right and may involve a great deal of modeling and 
data mining. When such information is not available, a 
simple choice, emulating the gravity model [•>], could be 
Vij = Avij , where 



J. 



(4) 



for i ^ j, and vu = 0. The idea here is that the at- 
tractiveness of a place, j, as a trip destination (from i), 
is proportional to its population, Cj, but inversely pro- 
portional to the travel distance from i to j, dij. The 
exponent 7 could be set equal to 2, acknowledging Klein- 
berg's condition for maximum navigability of the net- 
work through a decentralized algorithm [4], or to 7 = 1, 
so as to maximize the information entropy and the ease 
of collecting information in the social net ["1], but other 
choices might be equally reasonable. At this naive level, 
the same assumptions might be made for relocations, i.e.. 



The SMM can be simulated as a Markov process in a 
computer, using the rates of Table I, or it can be analyzed 
through a master equation, incorporating the very same 
rates. At the mean-field level, taking populations as con- 
tinuous variables and neglecting fiuctuations in number 



transition 



rate 



description 



{i, i) {j, i) 

(ij) {k,j) 
{i,k) 



When at home 



When away from home 

r 
r 

1 - r 



travel 
relocate 

return home 
continue trip 
relocate (optional) 



TABLE I: Transition rates for the SSM. In all cases, i / j / fc. 



space, the process is described by the rate equations: 



(5a) 
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(5b) 



optional 



We now turn to a particular implementation of the 
SMM that highlights the dramatic implications of its two 
modes of travel. Consider a basic susceptible-infected- 
recovered (SIR) model with births and deaths occurring 
at equal rate, /x, such that the total population, A'', re- 
mains fixed (on average). If all births are assumed to be 
susceptibles, then the governing equations are 



S^^fiN - ^iS~ (3SI, 



(6) 



along with the algebraic constraint that S + I + R = N. 
Here a describes the disease recovery rate and /3 is a 
transmission rate parameter. This model provides a rea- 
sonable representation for a disease such as the measles 
[()]. Standard analysis of (6) shows that if 
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> 1, 
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the disease is endemic and persists at some constant level, 
but if i?o < 1 the disease-free state is globally attracting. 
Assuming a specific disease and social setting, we may 
view a, /3, and /x as fixed and observe that whether the 
disease is endemic or dies out is determined by the pop- 
ulation size, N. 

Using this basic mean-field as a starting point, a num- 
ber of researchers have shown that a reasonable inter- 
pretation for many diseases is that large population cen- 
ters may operate at endemic levels, while outlying rural 
areas (where the disease would normally die out) have 
sustained levels or recurrent epidemics due to transport 
from the population centers [7] . Within that general con- 
struct, we examine a simple situation where we model 
a two-site system: site a, a city, and site b, a village, 
where (for specificity) we take the city population, Na, 
to be a factor of ten larger than the village population, 
Nb. Choosing appropriate parameters, we may consider 
a disease such that if two populations remain separate, 
the disease is endemic in the city and dies out in the vil- 
lage. But if we allow a small amount of social transport 
between locations, then the disease becomes endemic in 
both locations. We use this simple setting to show the 
dramatically different effects of the two types of motion; 
relocations and trips. 

Define our state variables as Xij , the number of people 
in health condition X G {S, I, R} who are at location 
i and call j their home, where i,j G Thus, for 

example, Saa are the number of susceptible people who 
are at the city a and a is their home. Applying a mean- 
field version of the SMM to (6) yields 

Saa = t-tXa — l-tSaa PSaa {laa + lab) 

— UJabSaa ~ l^abSaa + ^baSbb H Sba, (7) 

T 

which we describe and justify as follows: 

• All births occur while the person is home (not on 
a trip). 

• A susceptible is equally likely to get the disease 
from any infected person who is in the same loca- 
tion. 

• Relocations for all health conditions occur at the 
same rate Wafc (or tdba) and, following the gravity 
motion model, we further assume NaUJab ~ Nb-uiba- 

• Trips, for all health conditions, occur at rate 
i^ab (t'ha), again satisfying the gravity constraint 

Naiyab = NbVba- 

• We assume that relocations occur only from home. 

• The last four terms of (7) describe motion. For in- 
stance, the very last term refers to people returning 
from a trip (whose average duration is r). 



• We assume that the motions imply a change of loca- 
tion or home, but not a change of health condition. 

The rate equations for the remainder of the state vari- 
ables arc straightforward and follow the same principles. 
Rather than listing the full complement, we give just two 
further examples: 

laa — I^Saa {laa ^ lab) C^Iaa ^laa 

1 

—UJahlaa — Vablaa + (^bahb H Iba, (8) 

r 

and 

Sba = —fJ-Sba — fiSba {ha + hb) + ^abSaa Sbai (9) 

T 

where (8) re-emphasizes that transport preserves the 
health status, while (9) highlights the rule that no re- 
locations take place while not at home (state Xi^j^i). 
Wc now fix our parameter set to 

Na = 10, TVfc = 1, = 0.05, a = 48, /? 40, r = 5/365, 

which, according to the mean-field description, and in the 
absence of motion, supports an endemic level in the city 
a, but not in the village b. Apart from the average stay- 
ing time on a trip, r (which is fixed at r = 5 days), the 
motion is governed by the yet unspecified rates uJab and 
Vab- These can be constrained so that the total level of 
transport, from trips and relocations combined, is fixed 
at some small constant value. The relative volume of 
trips versus relocations is then controlled by a single pa- 
rameter 



that we term the population vacation index: If motion 
includes only relocations then 1/ = 0, and V grows ever 
larger as the fraction of trips increases. Our general con- 
cept that T > 6 would indicate that normally V > 1, 
though one can certainly imagine cultures where this 
might not be the case. 

Our SIR model has a single stable fixed point, which 
can be solved numerically. In Fig. 1 we show the steady- 
state one obtains in this fashion for the number of in- 
fected people in the village, as a function of the vacation 
index V. For "normal" conditions, ^ > 1, increased trips 
(offset by decreased relocations) results in a substantial 
increase of the disease level. For V < 1 (but increasing 
V), there is a decrease in the total number of infected 
individuals in the village. This results from the fact that 
sick people, coming from the city to the village on a trip, 
are less effective at transmitting the disease than those 
traveling for relocation, due to their much shorter stay. 
Trips become far more dangerous to the villagers when 
they don't compete with relocations {V > 1). 

In Fig. 2 we show the total fraction of susceptibles in 
the city and in the village (relative to their total popula- 
tions), under the same scenario. In the city, the large 
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FIG. 1: Epidemic levels in the village as a function of the 
vacation index V . Shown are Ihh and /ta, as well as the total 

lb = Ihh + lha- 




FIG. 2: Fraction of susceptibles, based at home. Observe that 
for low vacation index, the populations are essentially mixed, 
but when transport is dominated by trips, the locations are 
strongly heterogeneous, with the village having significantly 
greater risk of experiencing an epidemic. 



population acts effectively as a buffer, resulting in a 
nearly constant fraction of susceptibles even as the va- 
cation index changes from one extreme to another. The 
fraction of susceptibles in the village is nearly equal to 
that in the city, for lovif vacation index. However, as Y in- 
creases, there is a nearly four- fold increase in the fraction 
of susceptibles in the village, as trips become dominant. 

So far, we have analyzed the effect of T and B alone, 
ignoring the important role of r, the third characteristic 
time scale of motion in the SMM. In Fig. 3 we show the 
number of sick people in the village when the ratio of T 
to Q is held fixed, at V ~ 10, but the average duration 
of a trip, T, is varied from zero and up to 100 days. One 
can see that peak transmission occurs when the trip du- 
ration is roughly commensurate with the average length 
of the infection period, r « 1/a = 5 days. Shorter trips 
mean that infected individuals from the city return back 
home before fully spreading the disease, while longer trips 
allow villagers who travel to the city, and get infected 
there, to recover before returning home, thus reducing 
the level of disease in the village. Notice that this inter- 
esting "resonance" effect occurs only for trips, and could 
not happen when relocations predominate (as confirmed 
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FIG. 3: Change in village infection population as a function of 
trip duration. Fixing V — 10, we maintain transport volume 
while varying the trip duration, r. The resonance peak in the 
endemic level at r ~ 1/q occurs only for a sizable fraction of 
trips (see text). 



by additional numerical analyses). 

In summary, we have introduced a basic model for hu- 
man mobility that incorporates two very distinct modes 
of motion; trips and relocations, encompassing three 
widely different time-scales (two for trips and one for 
relocations). All of these time scales are potentially im- 
portant, as different phenomena might evolve at a charac- 
teristic time matching any of the three. Likewise, the fact 
that it is the very same individual that returns from a trip 
(as opposed to a randomly selected one) has important 
ramifications for the spread of disease. These hunches 
were demonstrated by a simple epidemic spreading SIR 
model evolving on just two locations, a "city" and a "vil- 
lage," showing that the endemic level found in the village, 
of a disease such as the measles, depends strongly on the 
mode of travel (whether trips or relocations). In addi- 
tion, a resonance effect, where the endemic level max- 
imizes, under unfavorable conditions, occurs only when 
trips predominate. Thus, making the distinction between 
the two modes of travel leads to very different numeri- 
cal results and reveals new qualitative phenomena. For a 
fixed level of "traffic" between location, we observe that 
the trips provide an increase in the disease coupling (com- 
pared to the relocation-only pattern) . It is possible that 
this effect may contribute the required coupling param- 
eter adjustment of a factor of 1.3 (as discussed in [:]]) to 
match England/ Wales measles data to a gravity trans- 
port model. 

Though couched in the specific language of human pop- 
ulations (and spread of disease, in our example), we an- 
ticipate that similar considerations would also apply else- 
where: animals, for example, might forage and return to 
their home base (trips), or leave the pack to form a new 
clan (relocation), on widely different time scales, and this 
might have implications to various aspects of their ecol- 
ogy, geographical genetic spread, etc. 

Our model, while hopelessly naive, incorporates the 
crucial differences between the two modes of travel — a 
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distinction that to date has gone largely ignored. Clearly, 
the model could be refined, by tweaking its rates, incor- 
porating other modes of travel, and so on, and such future 
revisions will result in ever greater predictive power. 
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